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Abstract 

We present a simple mechanical model for dynamic wetting phenomena. 
Metallic balls spread along a periodically corrugated surface simulating 
molecules of liquid advancing along a solid substrate. A vertical stack of balls 
mimics a liquid droplet. Stochastic motion of the balls, driven by mechan- 
ical vibration of the corrugated surface, induces diffusional motion. Simple 
theoretical estimates are introduced and agree with the results of the ana- 
log experiments, with numerical simulation, and with experimental data for 
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microscopic spreading dynamics. 



Typeset using REVT^ 



I. INTRODUCTION 



Because of its practical applications in areas such as coating, lubrication, adhesion, 
etc., the old field of wetting phenomena has recently attracted renewed interest Micro 
droplets which spontaneously spread along a solid surface have a time dependent shape 
which results from the balance between liquid - solid interactions and friction processes 
1^. Thus wetting mechanisms are intimately connected with friction on a molecular level 
|7| JT(]|JTT[| which is also important in the understanding of such practical problems as friction 
between two solid surfaces separated by a thin liquid layer |^ and dynamics of long polymer 
chains in random media [[T2 |. 



A salient feature of macroscopic spreading is that it is often preceded by a microscopically 
thin film - precursor. The precursor film thickness may vary from molecular size (one or 



sometimes several monolayers) to a few hundreds of angstroms. |J3[. For nonvolatile liquids. 



well below their critical temperature, thickness profiles with distinct successive molecular 



layers (terraces) have been observed |jT^JT^. EUipsometric measurements, carried out on 
different substrates and also for various kinds of simple liquids, as well as polymeric and 
surfactant melts, have reached a surprising conclusion: the linear size R of the precursor 



obeys a universal law p^ , p!6 | 



RocVi, (1) 

t being the time. The same law holds also for capillary rise, in which a vertical wall is put 
into a contact with a bath of liquid. Here a film of microscopic thickness grows from the 
macroscopic liquid meniscus and creeps upwards along the wall. In this case, the height of 
the film obeys the \/t-\aw within an extended time domain until it gets truncated,. 



at very high altitudes, by gravity. A diffusion-like coefficient Di can be formally defined as 
a prefactor in the Eq. ([l|), which is found to scale as the inverse of the bulk viscosity. Such 
a formal fitting does not, however, immediately imply an understanding or even a model of 
the microscopic dynamics responsible for the simple power law M. 
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Theoretical understanding of precursor dynamics has followed from two major concep- 
tualizations. The hydrodynamic approach (HA) |jT9[, and the Solid on Solid Model (SOSM) 
pO| which is based on Langevin equations for layers in the drop. Both approaches qualita- 
tively describe the formation of layered structures. HA correctly describes long term kinetics 
of terraced spreading in 2D systems with cylindrical symmetry, however the SOSM predicts 
R (X t rather theni/t. In both models the layers are considered as being incompressible con- 
tinua neither assumes nor implies a microscopic model of the dynamics. An analogy with 
the analysis of macroscopically thin layers [§, and Ising - like models @], as well as the 
diffusion like structure of Eq. (1), suggests that diffusion inside the precursor layer plays an 
essential role for spreading dynamics. However consideration of the motion of precursor edge 
(PE) as a simple biased random walk is incorrect 0. In the presence of an external force 
(capillary force which pulls the precursor out of the drop 0), the mean displacement of a 
biased random walker (at the PE) is proportional to time in violation of Eq. (1), as is the 



prediction for the dynamics of the first layer in the SOSM ||2^). In the absence of external 
forces the square root of mean squared displacement would exhibit the behavior of eq. (p. 
However, the mean PE displacement would be zero contrary to wetting experiments which 
indicate ^ a continuous directed displacement of the PE with relatively small fluctuations. 
This differs strongly from the fluctuation induced motion of a non-biased random walker, 
where the fluctuations are of the same order as the typical displacement. 



In [0,0 it was shown that the mean displacement of a random walker which is biased 
by a uniform external force and additionally experiences excluded volume repulsion exerted 
by an ensemble of other, non biased, diffusing hard core particles grows in proportion to 
v^, instead of the linear in time growth expected for similar systems without the hard core 
repulsion. We also have shown in this model that excluded volume effects imply an effec- 
tive frictional force imposed on the motion of an individual particle in the hydrodynamics 
description. 

The goal of the present work is to investigate a simple macroscopic spreading process with 
excluded volume. We consider the spreading of metallic balls, with and without magnetic 



interactions between the balls, and find that the spreading of balls in ID belongs to the 
same universality class as spreading of microscopic particles of liquids, i.e. is governed by 
eq. (|I|). We present experimental results for a the spreading of macroscopic balls, which 
emerge from a reservoir (vertical stack) , and spread onto a horizontal vibrating rack, see 
Fig. |I]. The surface of the rack is uniformly corrugated to prevent rolling or ballistic motion, 
of the balls. Constraints enforce zero or unit occupancy of corrugated sites, creating a lattice 
system. Weak random driving leads to lattice gas type behavior of balls - random jumps of 
length +/-1 constrained by the excluded volume effect. We find that this analog model yields 
the v^-law for the total number of balls which emerge from the reservoir, for the average 
displacement of balls from the reservoir origin, and for the displacement of the rightmost 
ball. We also present analytical estimates and numerical simulations, which are in good 
agreement with experimental data. The results suggest that, as far as time dependence is 
concerned, the \/i-law is essentially independent of the nature of the interactions between 
the substrate and the spreading substance, long range interactions between the particles 
themselves, as well as of the geometry and of the size of spreading particles. The prefactors 
in the square root law, of course, do depend on the system's parameters and all details of 
the microscopic interactions. 



The organization of the paper is as follows: In the Sections |ll7\| and |IIB| we describe 



the experimental set-up and numerical algorithm, in the Section 111 we present a simple 



analytical results for our model: basic equations - in the Section pi7\| , results for the flux of 
balls from the reservoir and averaged displacement of spreading balls in the Section |III B| , 
results for the displacement of the rightmost ball, that determines the size of the spreading 
layer - in the Section [III C| , and separate results for the steady regime which starts after the 
firs ball falls out of the rack - in the Section |111 D| . In the Section jTV] we present experimental 



and numerical results, which are briefly summarized in the Section - |V|. 



II. EXPERIMENTAL SET-UP AND NUMERICAL PROCEDURE 
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A. Experiment 



A schematic diagram of the experimental set-up appears in Fig. [l|. A corrugated 
(notched) horizontal rack is confined inside a rectangular tube which prevents balls from 
passing one another or jumping off the track. A vertical stack is placed at the left end of 
the rack (at the origin) and metallic balls (balls) are fed through this stack maintaining 
unit concentration of 'particles' at the origin. The balls are allowed to move to the right 
of the origin only. The entire system is driven with motors placed on each end of the rack. 
The fiywheel on each motor is eccentric to provide "chaotic" oscillations which are exper- 
imentally shown (see below) to give rise to diffusional motion of the balls. A particular 
number of balls, n (n = 4 or 8), is placed in the vertical stack and then the oscillations are 
started. The time for the balls to leave the stack is measured and recorded along with the 
displacement of the horizontal balls. After this, n more balls are added to the stack and the 
procedure is repeated with the first n balls left at their respective place on the horizontal 
rack. The process is repeated until the rightmost ball reaches the right end of the rack. 
Results presented are averaged with respect to four independent trials. In a separate set 
of experiments the number of balls which emerge from the vertical stack as a function of 
time after the rightmost ball reached the right end of the rack and escapes from the rack 
was measured. This is the spreading rate for a "full" horizontal rack. Two kinds of balls we 
used: magnetic and non magnetic. 

The experiment was designed to mimic the spreading of liquids: the vertical stack imi- 
tates a liquid drop, which acts as a reservoir providing particles for the precursor and keeps 
constant concentration at the drop - precursor boundary, the horizontal rack mimics a solid 
surface. Since the gravitation energy of balls on the rack is proportional to the height, the 
shape of the surface of the rack emulates the profile of the potential of liquid - substrate 
interactions with local potential minima, which prevent particles from long range (on length 
scales larger then the mean distance between two neighboring particles) ballistic or rolling 
motion. Oscillations of the horizontal rack mimic thermal excitations of the surface of a solid 



6 



substrate and generate random jumps of balls along the rack. The number n determines the 
gravitation induced pressure in the vertical stack. Two sets of experiments - with n = 4 and 
n = 8 were carried out in order to make sure that in given range of n's values the pressure 
is not important, see section 0. The total number of balls on the horizontal stack, M{t), 
is an analog of the precursor mass and the displacement of the right most ball is an analog 
of the precursor radius, R{t). The magnetic mutual ball to ball attraction serves to mimic 
the particle - particle attraction in a precursor film and in a liquid drop. 

To test the assumption that the eccentric oscillations generate diffusive motion of a 
single ball. Fig. ^ shows the mean squared displacement and the squared mean displacement 
of a single ball initially placed on the horizontal rack half way between the origin and 
right-most end. The observed linear dependence of the mean squared displacement on 
time agrees with well known result of conventional single particle random walk theory. 
The observed RMS dispersion from the mean is also in approximate agrement with the 
prediction of simple estimates presented in the Appendix. However, the absolute values of 
dispersion is approximately 1.5 times greater then the theoretical value. These experimental 
result indicates that the jumps of a single ball are, to a good approximation, random and 
independent events caused by the vibrational driving of the system. 

B. Numerical simulations 

In numerical simulations, we modeled random jumps of particles on a ID lattice with 
unit steps. For each time, each particle except the particle in the site number chooses 
randomly a direction: to the right or to the left with equal probabilities. If the corresponding 
neighboring site is occupied, the move is rejected and the particle stays at the original 
position; if the site is vacant, the particle jumps to it. The origin (site 0) is always occupied; 
a new particles is automatically added to the site number 1 when it becomes vacant. To 
simulate an analog of magnetic interaction in the simplest manner, the jump rates of particles 
which had a nearest neighbor from one side to the vacant site on the other side was reduced by 
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the factor 2(1 —p) where p determines the "strength of the interaction", p = 0.5 corresponds 
to the absence of the interactions and p = 1 corresponds to infinitely strong attraction. At 
an initial time, two particles are placed into the system: the first particle - at i = and the 
second particle - at z = 1. The mean displacement of balls was measured as was the number 
of balls as function of time up to 100 active balls on the rack. Results were averaged with 
respect to 40 independent trials for each p = 0.5, 0.6, 0.7, 0.8. 



III. SIMPLE ANALYTICAL ESTIMATES 

A. DifFusion equation and boundary conditions 

In order to obtain a simple analytical estimate for the dynamics of the processes we 
neglect mutual magnetic attraction (this is the same as neglecting surface tension for mi- 
croscopic spreading phenomena), and decouple many particle probabilities. This leads to 
following equations for mean concentrations in ID lattice with unit step (horizontal rack): 

^^ = ^{C,(^ + l,t)[l-C,{^,t)] 
+C,{t - 1, t) [1 - C,{t, t)] - C,{t, t) [1 -C,{t- 1, t)] , (2) 
-C,{t,t)[l-C,{t + l,t)]} 

where to is frequency of jumps, which is connected to the diffusion coefficient, D = I is 
the jump length, Cv{i,t) - concentration of vacancies (empty slots on the horizontal rack) 
at the slot number i at time t, Cy{i,t) = 1 — Cb{i,t), Cb{i,t) is the probability to find a 
ball (ball) at slot number i at time t. The first term in the right hand side of the Eq. (^ 
determines the rate of jumps of vacancies from the site number i + 1 to the site number i. 
It is proportional to the concentration of vacancies at the site number i + 1 multiplied by 
the concentrations of particles at the site number i since only vacancy - particle exchanges 
are allowed. The other terms determine the rate of jumps i — 1 i, and z — > i — 1 and 
i — * i — 1 in a similar way. 
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In the approximate Eq. (j^) , nonlinear terms cancel since the forbidden particle - particle 
and vacancy - vacancy exchanges do not alter the local mean concentrations. When i and t 
are large, z >> 1 and ut » 1, it leads to the diffusion equation: 

^ = |aC.,(M). (3) 

where A is Laplace operator with respect to the variable i. At t = and i > 0, there 
is no "liquid" on the solid surface (no balls on the horizontal rack). Therefore, the initial 
concentration of vacancies is equal to unity: 

a(z,t)|,^o = l (4) 

- all sites are empty. When i is large the concentration of vacancies is equal to unity because 
the balls have not had enough time to reach this area, 

lim C^(i,t) = 1. (5) 

There are no vacancies at the boundary i = 0, because the reservoir of balls places a ball in 
the vacancy instantly, 

a(2,t)Uo = 0. (6) 

Note that in spite of the fact that the Eq. apparently does not reflect excluded volume 
interaction, the boundary condition Eq. states that the ball from the vertical stack can 
move down to the horizontal rack if and only if the vacancy comes to the point i = 0. When 
the ball falls down from the stack it eliminates the vacancy at t = 0. The relaxation time for 
the vacancy concentration in the horizontal stack is determined by the diffusion of a vacancy 
through the array of balls in the horizontal stack - from the right end to the left end. Since 
the number of balls and the length of the array grow, the relaxation time increases which 
slow down the dynamics and leads to the dependence presented by Eq. (|1|) instead of linear 
growth. 
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B. Flux from reservoir, number of balls, and averaged displacement 



The solution of the Eq. (H) in ID with the boundary conditions eq. (M and eq. 



C,{t,t) = l-erfci—^) 



2VDt 

leads to the following result for the flux, P{t), at i = 



(7) 



P{t) 



2 di 

The number of particles on the rack is equals to 



i=0 



27rt' 



(8) 



M{t) = J P{T)di 



l2ujt 



(9) 



Writing down the equation similar to Eq. for particles concentration, multiplying both 
sides of this equation by i and integrating from to infinity, and taking into account the 
boundary conditions Eq. and Eq. (H), we obtain, for the mean total displacement of all 



balls in the horizontal rack, K{t) = J iCb{i,t)di, measured in the jump length units / 





dK{t) 
dt 



1 

-OJ. 

2 



(10) 



Integrating Eq. ([T0|) and taking into account the normalization 

oo 

M{t) = I Cb{t,t)di, 



we obtain the averaged over all balls displacement for the balls in the horizontal rack, r{t) 

Kit) 1 



r{t) 



-V2TCUjt. 



(11) 



M(t) 4 

The straightforward extension of the results, Eq. and Eq. (|TT1) for the corresponding 
2D system in large t limit leads to 



Pit)\ 



2nau! dCv{i, t) 



2D 



dr 



oc 



i=0 



21n(t?) 



(12) 



and 



M(t)|2^cxr(t)| 



TfUjt 



OC 



^^-21n(f)- 



(13) 
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C. Displacement of the rightmost ball 



One can show that apart from logarithmic corrections, which can occur in long time 
regime of spreading, the displacement of the rightmost ball on the horizontal rack (the length 
of the precursor- R{t)) is proportional to the average displacement, R{t) oc r(t), which leads 
to the Eq. (|I|). Here we present a simple estimate for the R{t). The distribution of the 
displacement of the right most particle, for the ensemble of spreading hard core particles 
can be bounded by the maximum displacement of independent diffusing particles. The later 
displacement has the following distribution 

P{R) = c(i?,t)e-^ri'^(^-'=("'*))''^ (14) 

where c{R,t) = 1 — C^(i,t), C^{i,t) being determined from the Eq. (|^. The first multiplier 
in the right hand side of the Eq. (|ri|) is the probability to have a particular displacement 
R for a ball, the second term is a limiting continuum form for a probability to have smaller 
displacement for other balls. Averaging R with respect to the distribution Eq. (|1^) by means 
of steepest decent method in large t - limit we obtain 

V vra^ 
D. Steady regime 

The regime of spreading changes after the first ball falls off of the horizontal rack. The 
length of the array of balls does not grow any more and the concentration of balls is de- 
termined by the steady state solution of the Eq. (H) with a new boundary condition at the 
right side of the horizontal rack 

C^.(^,'^)I.=L = 1, (15) 

where L is the total number of slots on the horizontal rack. This boundary condition states 
that there are no balls out of the rack with the coordinates i > K. The solution of the 
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steady state Eq. leads to the constant flux of vacancies to the origin, which in turn leads 
to 

Mi(t) = ^(t-to), (16) 

where Miis the number of balls which fall out of the rack and to is the time when the first 
ball falls. 

IV. EXPERIMENTAL AND NUMERICAL RESULTS 

The Fig. ^ and Fig. ^ show the experimental results for the number of balls in the 
horizontal rack and the averaged displacement of the right most ball versus square root of 
time for magnetic and non magnetic balls for n = 4, 8 and t < to, 

Fig. ^ shows the time for a given number of balls to enter the horizontal rack from the 
stack . In this case the experiment was not interrupted at time t = to. 

The Fig. ^ presents the dependencies of the mean squared displacement along with the 
squared mean displacement of the rightmost ball in the horizontal rack. The figures Fig. 
Fig. 1^, and Fig. ^ show the dependencies of the mean number of balls in the horizontal rack, 
mean displacement of balls on the horizontal rack, and the mean displacement of the right 
most ball versus the square root of time obtained in the numerical simulation for different 
values of the "interaction parameter" p. 

The pronounced straight lines which are presented in Fig. ^ Fig. |, Fig. Fig. |^, and 
Fig. H are in good qualitative agreement with Eq. (pJ]) , Eq. @ and Eq. ([l|). However, 
the fluctuations in the dependence of the mean displacement on square root of time are 
larger then fluctuations for the dependence of the total number of balls. The results for 
n = 4 and n = 8 are essentially similar. It shows that the gravity - induced pressure in the 
reservoir does not play significant role, which corresponds to the spreading of small droplets 
or vertical creep with relatively small precursor length where the gravitational forces are not 
important. In a test experiment with larger number of balls, n > 40, the large gravitational 
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force pushed all balls from the vertical stack after beginning of vibrations, overcoming the 
potential barriers produced by slots on the horizontal rack, i.e. large gravitational forces 
changed the nature of spreading. 

The theoretical values for the slopes of the linear dependencies of the number of balls 
and of the mean displacement of balls in the precursor which are determined by the Eq. 
and Eq. (|ll]) for = 1: ^Jl ^ 0.798 and ^ 0.6266, are in excellent agreement with 

the corresponding values determined from the numerical experiment : 0.804 and 0.625. 

Magnetic interaction does not change the shape of the dependencies but decreases the 
numerical prefactors and increases fluctuations. Introduction the effective interaction in 
numerical simulations leads to similar effects. The numerical results for the dependencies of 
the mean number of balls on the horizontal rack and of the mean displacement versus square 
root of time, presented in Fig. ^ and Fig. |] exhibit a slowing down when the "interaction 
parameter" increases from 0.5 to 0.8. The fluctuations in the mean squared displacement 
also increase and crossover period, which can be seen for small times in Fig. § becomes 
larger. 

It is also instructive to compare the results presented in Fig. ^ for the moments of the 
displacement of the single ball on the horizontal rack with the results for the moment of 
the displacement of the rightmost ball of the array spreading along the rack, presented in 
Fig. 1^. For a single ball the mean squared displacement is proportional to time while the 
squared mean displacement is much smaller and irregular. In the limit of large number of 
trials, it should tend to zero as l/V^On the contrary, the experimentally measured 
mean squared displacement and squared mean displacement for the rightmost ball in the 
array are equal within the error of the experiment. This shows that, in spite of apparent 
scaling, the similarities of these two processes are of essentially different physical origin and 
behavior. The motion of a single ball is a fluctuation induced process with zero mean, while 
the spreading is a driven diffusive process with small fluctuations. 

The initial parabolic dependence presented in Fig. 5 for t < to agrees with the Eq. 
while the linear dependence obtained for t > confirms Eq. (|T6|). The numerical prefactors 
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for these time dependencies are not the subject of our present work; thus actual values of 
the vibrational frequencies and length scales are not important for our results. However, 
we can suggest some rough estimates, which produce reasonable numerical values. The 
effective experimental frequency of jumps uj ~ 32sec~^ was roughly estimated from the data 
presented in Fog. 5 for the time t > to (the liner regime) with the approximate values for the 
length of the rack ^ 60cm and / ~ 1cm. It is in good agreement with the values u ~ 36sec~^ 
which was obtained by independent estimate from the region t < to my by means of Eq. @. 
The slope for the experimental dependence of mean displacement versus square root of time 
is also in good agreement with the theoretical estimate. 



V. SUMMARY 



We have calculated analytically, simulated numerically, and measured experimentally 
the mass and the size of a wetting layer for a simple mechanical system - metallic balls, 
which spread along a periodically corrugated surface, which mimics the wetting phenomena. 
The experimental results and numerical data are in excellent agreement with the theoretical 
predictions concerning the dependencies which determine the dynamics of spreading of a 
lattice gas (balls) in our mechanical model . These agree with Eq. (|l|) which empirically 
describes the dynamics of precursors in spreading of actual liquid drops. The mechanical 
model may thus reflect the important features of this phenomena which lead to the universal 
spreading law, eq. ([l|). Analysis of the experimental results, simple analytical estimates 
and numerical simulations provides an additional argument that the precursor kinetics is 
controlled by the diffusion of vacancies from the precursor's boundary to the liquid drop 
H^,|7|p^. A theoretical model which directly explores the consequences of these ideas for 



a model of microscopic wetting is presented in a subsequent paper [Q, which substantially 
extends the ideas of [E^ . 
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APPENDIX A: DISPERSION OF MEAN SQUARED DISPLACEMENT 

Let us consider the dispersion for the average over K realizations squared displacement 
of a random walks containing A'" steps. What we measure is: 

1. mean squared displacement for individual random walk 



where sj is the step number j, Ri{N) is the averaged squared displacement of the 
random walk number i. 

2. squared displacement for individual random walk realization 



4. square root from averaged deviation of mean squared displacement from it's averaged 
value 





3. averaged with respect to all realizations squared displacement 




K 



K 



( 



K 




{N) 



K 



i=l 




V 



K 



Denote Rf 
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All Ri{N) are independent random variables with zero mean and second moment 
equal to Na"^, the gaussian random variables Rf{N) have mean equal to iVcr^, and sec- 
ond central moment equal to 3A^V^ For K = 30,^30 oc 5.4772, 1/^30 = .18257 and 
A = ^Js/V^Na'^ = .74008 A^cr^; 
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FIGURES 



FIG. 1. Sketch of experimental set-up. 
FIG. 2. Mean squared - diamonds, and squared mean - triangles, diplacement of a single ball 

vs time, averaged over 30 trials. 

FIG. 3. The experimental dependence of mean number of balls on square root of time. Triangles 

- non magnetic balls, n = 8; squares - non magnetic balls, n = 4; diamonds - magnetic balls, solid 

lines - linear regression. 

FIG. 4. Mean displacement of the rightmost ball vs square root of time. Triangles - experimen- 
tal data for non magnetic balls, n = 8; squares - experimental data for non magnetic balls, n = 4; 

diamonds - experimental data for magnetic balls, solid lines - linear regression. 

FIG. 5. Time vs number of balls which emerged from vertical stack. Squares experimental data 

for non magnetic balls. Solid line - theoretical dependence predicted by Eq. (8) for t <to, dashed 

line - theoretical curve, predicted by Eq. (12). 

FIG. 6. Mean square displacement - squares, and square of mean displacement - triangles, vs 

time for n = 8 and non magnetic balls. 

FIG. 7. Numerical data: mean number of balls on the horizontal rack, p grows from 0.5 to 0.8. 
FIG. 8. Numerical data: mean displacement averaged for all balls on the horizontal stack, p 

grows from 0.5 to 0.8. 

FIG. 9. Numerical data: mean displacement of the rightmost ball along the horizontal stack, p 
grows from 0.5 to 0.8. 
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